Lecture 14

class: center, middle, inverse, title-slide

# Lecture 14
## Models for Factorial Designs
### Psych 10 C
### University of California, Irvine
### 04/29/2022

---

## Review

- This week we have been going over factorial designs and the models that we 
use to analyze the results of such experiments.

- The first model we talked about was the Null model, which assumes that the 
expected value of the dependent variable is the same regardless of values that
a factor can take.

- We also talked about the Main effects models, which assume that only one 
factor has an effect on the expected value of the dependent variable. These 
models are separated by factor, which means that there are only differences 
between the groups that respond to the levels of one factor regardless of the 
values of the other.

- Additive models assumed that 2 or more factors could have an effect on a 
participants' response, however, that effect was considered to be independent 
and therefore was equal to the sum of the main effects of each factor `$(\alpha$` 
and `$\beta)$`.

- Finally, we talked about the Full model, which assumes that the expected 
value of the dependent variable was different for each group (combination of
factor levels) in an experiment. Furthermore, this expectation can't be 
expressed as the sum of the independent main effects of each factor.
 
---

## Review

- The prediction of the Null model in factorial designs is known as the 
grand mean, and it's the average of the responses of all participants regardless
of the level of each factor that they where exposed to. `$$\hat{\mu}$$`

- The predictions of the Main effects models where equal to the sum of the 
average response across all participants and the "effect" of the factor that 
defines each model.$$\hat{\mu} + \hat{\alpha}_j$$ or `$$\hat{\mu} + \hat{\beta}_k$$`

- For factor 1 or factor 2 respectively.

---

## Review

- We mentioned that using the estimates of the main effects `$\hat{\alpha}_j$` and
`$\hat{\beta}_k$` we can derive the predictions of the additive model as:
`$$\hat{\mu} + \hat{\alpha}_j + \hat{\beta}_k$$`

- Finally, we mention that the predictions of the full model where equal to the 
average of our dependent variable for each combination of the levels of our 
factor: `$$\hat{\mu}_{jk} = \frac{1}{n_{jk}} \sum_i y_{ijk}$$`

- Where `$n_{jk}$` represents the number of participants in the group that 
responded to the combination of the *j-th* and *k-th* values of factors 1 and 
2 respectively.

---

## Cell means

- The equations that we have seen in previous classes will always allow us to 
obtain the values of the parameters in each model. However, when we have a 
between subjects factorial design we can use the cell means method in order
to obtain those values.

- This method requires us to first calculate the average of each combination
of the levels of our factors.

- When we have a small number of participants on each group this is easy to do 
by hand. We will try this today.

---

## Cell means mehtod

- First we need to introduce the problem that we want to solve.

- We want to know if the cohort that a student belongs to or whether or not they
took a statistics class during their first year has an effect on the anxiety 
levels of students at the end of their first year.

- In order to study whether our factors (cohort and statistics class) have an 
effect on the anxiety levels of students we asked students at a university to
respond to an anxiety scale.

- Then we recorded whether they started college in 2019 or if they started in 
2020 and whether they took a statistics class during the current year or not.

- We collected information of a total of 16 participants for this study, 4 of 
whom started college in 2019 and took a statistics class, 4 started college in 
2019 and didn't take a statistics class, 4 started college in 2020 and 
took a statistics class and the remaining 4 started in 2020 and didn't take a 
statistics class during the last year.

---

## Anxiety data:

<div id="htmlwidget-753ccab0020033b9c5ba" style="width:100%;height:auto;" class="datatables html-widget"></div>
<script type="application/json" data-for="htmlwidget-753ccab0020033b9c5ba">{"x":{"filter":"none","vertical":false,"fillContainer":false,"data":[["1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16"],["1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16"],[9,8,8,9,10,12,12,13,11,13,12,13,14,15,15,16],["2019","2019","2019","2019","2019","2019","2019","2019","2020","2020","2020","2020","2020","2020","2020","2020"],["other","other","other","other","stats","stats","stats","stats","other","other","other","other","stats","stats","stats","stats"]],"container":"<table class=\"display\">\n  <thead>\n    <tr>\n      <th> <\/th>\n      <th>id<\/th>\n      <th>anxiety<\/th>\n      <th>cohort<\/th>\n      <th>class<\/th>\n    <\/tr>\n  <\/thead>\n<\/table>","options":{"pageLength":8,"columnDefs":[{"className":"dt-right","targets":2},{"orderable":false,"targets":0}],"order":[],"autoWidth":false,"orderClasses":false,"lengthMenu":[8,10,25,50,100]}},"evals":[],"jsHooks":[]}</script>

---

## Anxiety example

- Using the data from slide 7, obtain the cell means for the anxiety problem.
Remember that the cell means are the predictions of the Full model, in 
other words, they are equal to the average of each combination of the factor
levels.

|          | Statistics            | Other                 | Mean                 |
|----------|:---------------------:|:---------------------:|:--------------------:|
| **2019** | `$\hat{\mu}_{11}$`      | `$\hat{\mu}_{12}$`      |  |
| **2020** | `$\hat{\mu}_{21}$`      | `$\hat{\mu}_{22}$`      |  |
| **Mean** |  |     |  |

- Fill in the following values with the correct average.
<br>

.pull-left[
|      | Statistics  |
|------|:-----------:|
| 2019 | .can-edit[] |
| 2020 | .can-edit[] |
]

.pull-right[
|      | Other       |
|------|:-----------:|
| 2019 | .can-edit[] |
| 2020 | .can-edit[] |
]

---

## Anxiety example

- After calculating the cell means, we can now calculate the mean by factor, for
example, the mean of students in the 2019 cohort would be equal to the average 
of the values on the first row, and for the 2020 cohort it is equal to the values
on the second row.

|          | Statistics            | Other                  | Mean                 |
|----------|:---------------------:|:----------------------:|:--------------------:|
| **2019** | `$\hat{\mu}_{11}$`      | `$\hat{\mu}_{12}$`       | `$\hat{\mu}_{1\cdot}$` |
| **2020** | `$\hat{\mu}_{21}$`      | `$\hat{\mu}_{22}$`       | `$\hat{\mu}_{2\cdot}$` |
| **Mean** | `$\hat{\mu}_{\cdot 1}$`  | `$\hat{\mu}_{\cdot 2}$` |  |

- Fill in the following values with the correct average.
<br>

.pull-left[
|      |  Mean       |
|------|:-----------:|
| 2019 | .can-edit[] |
| 2020 | .can-edit[] |
]

.pull-right[
|            | Mean        |
|------------|:-----------:|
| Statistics | .can-edit[] |
| Other      | .can-edit[] |
]

---

- Finally we can compute the grand mean by taking the average of the first 
cell means (the predictions of the full model).

- The grand mean was equal to :

- `$\hat{\mu}=$` .can-edit[]
  
---

## Anxiety example

- With the values that we have obtained we can now get the main effects 
`$\hat{\alpha}_j$` and `$\hat{\beta}_k$`. That way we will be able to get the
predictions of the additive model.

- To get the value of `$\hat{\alpha}_1$` remember that we need to get the 
difference between the the grand mean and the average of the group.

- Fill in the following table with the values of the main effects:

|                  |  Main effect |
|------------------|:------------:|
| `$\hat{\alpha}_1$` | .can-edit[]  |
| `$\hat{\alpha}_2$` | .can-edit[]  |
| `$\hat{\beta}_1$`  | .can-edit[]  |
| `$\hat{\beta}_2$`  | .can-edit[]  |

---

## Predictions by model

- What are the predicted anxiety levels of each model?

- Null

- Main effects Cohort

---

## Predictions by model

- Main effects Class

- Additive

---

## Predictions by model

- Full